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Abstract 

A systematic method to obtain strong coupling expansions for scattering quantities in Hamiltonian 
lattice held theories is presented. I develop the conceptual ideas by means of the Hamiltonian 
held theory analogue of the Ising model, in d space and one time dimension. The main result is a 
convergent series representation for the sacttering states and the transition matrix. To be explicit 
the special cases of d = 1 and d = 3 spatial dimensions are discussed in detail. I compute the next- 
to-leading order approximation for the phase shifts. The application of the method to investigate 
low-energy scattering phenomena in lattice gauge theory and QCD is proposed. 



1 Introduction 



The investigation of scattering phenomena in lattice field theories has attracted increasing interest 
in recent years. Since the work of Liischer [1], who established a close relation between the two- 
particle spectrum in a sufficiently large volume and the elastic scattering amplitude, several authors 
have succeeded to extract scattering data from numerical simulations of lattice models in 2,3 and 4 
space-time dimensions [2]-[5]. As further steps towards the determination of the low-energy scattering 
parameters of QCD there have recently been encouraging attempts to investigate resonance phenomena 
[6]-[8], and even some first full QCD calculations have been performed [9, 10]. 

This paper advocates a different method to study scattering in lattice field theories. In particular I 
present a systematic technique to obtain strong coupling expansions for the scattering matrix and the 
phase shifts in the elastic regime, applicable to a large class of lattice models including gauge theory 
and QCD. Such a proposal may look a bit old-fashioned. But if one remembers that strong coupling 
expansions have provided important qualitative insight into numerous issues of lattice field theories 
it is astonishing that only little effort has been made to apply this powerful analytical method for a 
computation of scattering quantities [11]— [13]. To close this gap towards a thorough understanding of 
scattering in lattice field theories is the intent my approach. 

Since scattering processes involve asymptotic two-particle states whose real time evolution is not 
directly accessible on a euclidean lattice I prefer the Hamiltonian formulation to develop the basic 
ideas. For many lattice field theories the Hamiltonian consists of two essentially different terms 



where k is inversely proportional to some coupling constant. The first part, IH , is a sum of terms 
which refer to a single lattice site, whereas IH x couples neighbouring sites. As long as the inverse 
coupling k is zero the system is purely static and it can usually be solved, if not analytically, then 
numerically. Turning on the parameter k is what makes the dynamics set in. Actually IH x creates 
and annihilates lattice excitations at different sites and thereby induces their propagation through 
the lattice. One may try to describe the scattering of those excitations by means of a perturbative 
expansion in the inverse coupling. 

As the test case I choose the Hamiltonian field theory analogue of the Ising model, in d space and 
one time dimension. This theory is obtained by considering the transfer matrix of the Ising model in 
the extremely anisotropic limit of the exchange couplings [14]. It is complicated enough to confront 
us with most of the conceptual difficulties of the problem. Still the model is sufficiently simple (no 
resonances, for example) such that the basic issues are not obscured by technicalities. 

The organization of this paper is as follows. To begin with I discuss the static limit of the (d+ 1)- 
dimensional Ising model. After briefly recalling a systematic perturbation theory for isolated levels 
this method is applied to the static vacuum, the one-particle and two-particle sector. The outcome 
are effective Hamiltonians which determine the dynamics of the system in these sectors. The rest of 
the paper is devoted to the diagonalization of the effective two-particle Hamiltonian. In section 3 I 
describe the general method to contruct convergent series representations for the scattering states and 
the transition matrix. The leading and next-to-leading order contributions are calculated explicitly. 
I continue by studying the simplest case, the (1 + l)-dimensional system in some detail. In section 
4 the next-to-leading order approximation for the scattering states is determined in closed analytical 
form such that the physical properties can be investigated. Life is substantially more complicated in 
the case of d = 3 spatial dimensions which is considered in section 5. At small momentum the con- 
tinuum rotational symmetry is restored on the lattice. Thereby the diagonalization of the scattering 
matrix simplifies considerably . By means of a transformation to better suited momentum variables 
the analysis can be extended beyond the small momentum regime. The next-to-leading order approx- 
imation for the scattering phases is computed. This paper ends with a few concluding remarks and 
is completed by three appendices, two discussing the physical properties of the one-particle and the 
scattering states and one discussing the calculation of the three-dimensional Green's function. 



IH = IH + kIH 
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2 Effective Hamiltonians 



The system studied here is the Ising model in (d+ 1) dimensions. In the static limit where the inverse 
coupling k is zero the spectrum consists of isolated eigenvalues. This is the prerequisite for a systematic 
perturbative expansion of the physical quantities by means of a very elegant method developed by 
Bloch [15]. The leading and next-to-leading order correction to the static vacuum, one-particle and 
two-particle sector are calculated. 

2.1 The ( d + l)-dimensional Ising model and its static limit 

We consider a d- dimensional hypercubic spatial lattice 

r = {x | x/a e TL\ -L/2 < x k < L/2} (2.2) 

with linear extent L and spacing a which is chosen equal to 1 for convenience. So L is an integer 
number and the physical quantities are dimensionless. Furthermore we assume periodic boundary 
conditions. At each site x of T we associate a copy of C 2 and call it C 2 . The unit vectors 

(2.3) 

are referred to as spin up and spin down states at x, respectively. For finitely many spins a Hilbert 
space 7i is defined as the tensor product of all C 2 . Let <7 1 ,<7 2 ,0'3 denote the familiar Pauli matrices 
and (7i(x) represent the operators on 7i that act by cr, on the x component of the tensor product and 
by the identity on the others. These operators commute at different sites which guarantees that the 
corresponding spins are independent observable quantities. 

The field theory version of the Ising model, in d space and one (continuous) time dimension, has 
the Hamiltonian [14] 

IH' = |{IH + ^IHJ , (2.4) 

where 

IH = ]T{1I - a 3 (x)} , (2.5) 

d 

IHi = -^^^(x^x + k) . (2.6) 

xer k=i 

The parameter g is a dimensionless coupling constant, k = 2/g 2 and k denotes the unit vector in the 
positive ^-direction. For the following it is convenient to define 

IH = - IH' = IHo + KlHi . (2.7) 
9 

As already mentioned in the introduction the Hamiltonian consists of two essentially different terms. 
The static part IH does not couple different lattice sites. On the contrary, a x is the operator that 
flips an up spin down and a down spin up. So the hopping terms (Ji(x)ai(x -\-k) flip the spins on 
neighbouring sites and thereby propagate signals through the lattice. 

The idea of this work is to investigate scattering phenomena by means of a perturbative expansion 
in the inverse coupling k. To begin with we have to determine the spectrum of the unperturbed static 
Hamiltonian IH . The k = ground state is determined through the condition 

IH |O)° = , (2.8) 
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i.e. the strong-coupling vacuum is the unique state having all spins up. It is normalized such that 
°<0|0>° = 1. We define 

\x)° = CTi(a;)|fi) (2.9) 
as the state where the spin is down at site x and all others are up. One finds 

\H \xf = 2\x)° , (2.10) 

hence (Ji(x) creates a one-particle state of energy 2. The successive action of a x leads to multi-particle 
configurations 

\x\...,x n )° = a^x 1 )..^^)^) (2.11) 
which are the higher excitations of the static system 

IHola; 1 ,...,^) = 2n\x\...,x n )° . (2.12) 

Since the spin-flip operators at different sites commute, any permutation of {x 1 , . . . , x"} in the above 
definition (2.11) represents the same ra-particle state. But it has to be emphazised that a x and 
a\ = Gi are no bosonic creation and annihilation operators in the usual sense. Rather they are 
unitary operators satisfying the constraint a\ = U which leads to (T 1 (a;)(T 1 (a;)|i7) = That is, two 

excitations cannot occupy the same lattice site and the multi-particle states are characterized by a set 
of mutually different vectors in 7L d . So the ra-particle eigenspaces of IH are given by 

E n = span {(a; 1 ,..., a;") | x 1 ^ ... ^ x n ,x k G T} . (2.13) 

They are well separated by an energy gap of 2 and highly degenerate. 

2.2 Bloch's perturbation theory 

We are concerned with the self- adjoint operator IH = IH + k Mi defined on the Hilbert space 7i. The 
static part IH has a purely discrete spectrum with eigenvalues 2ra. The projection operators on the 
corresponding eigenspaces 8 n are denoted by P„. 

If the perturbation k\Hi is turned on, one expects that each static eigenspace 8 n is deformed into a 
linear space 7i n C TC. In 7i n there is a complete (discrete or continuous) set of states \n,a) such that 

IH |ra, a) = E n (a)\n, a) , E n (a) = 2ra + e n (a) , e n (a) = 0(k) . (2-14) 

To obtain the eigenstates \n,a) and the eigenvalues e n (a) as an explicit power series of the inverse 
coupling Bloch proceeds in two steps. First the Hamiltonian IH is reduced to a hermitian operator 
acting on 8 n whose eigenvalues are equal to e n (a) to all orders of k. In a second step this effective 
Hamiltonian is diagonalized. Let |ra,a)° denote the eigenvectors corresponding to the energy values 
e n (a). Then the announced expansion of the eigenvectors of IH is obtained by means of a linear 
operator U n which maps |ra,a)° onto \n,a). 

Both and U n are explicitly constructed in Bloch's paper [15]. If the inverse coupling is suffi- 
ciently small they are well-defined power series of k. Each summand is a product of the perturbation 
IH ! and the projection operators P n and Q n = U — P n . In this work we only need the leading and 
next-to-leading order terms of the hermitean reduced Hamiltonian 

H; = «P„IH 1 P„ + ^PJH, Qn IH^ + 0(k 3 ) . (2.15) 

In — IMq 

and the linear operator 

U n = P„ + k Qn IH^ + C^k 2 ) . (2.16) 

In — Inn 
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Figure 1: Elementary graphs. 

2.3 The vacuum sector 

The vacuum eigenspace £ is the one- dimensional subspace of 7i with projection operator P = 
|fi) 00 (fi|. The action of the perturbation 

d 

IH i = -'J2'J2a 1 (x)a 1 (x + k) (2.17) 

k = l x 

on the strong-coupling vacuum creates two-particle excitations with flipped spins at neighbouring 
sites. From (2.15) it follows that the leading order contribution to the reduced Hamiltonian H' 
vanishes since states from different static eigenspaces are matched. In the next-to-leading order term 
Qo represents the sum over all intermediate states orthogonal to £ . But here only the two-particle 
states contribute for the same reason. Hence Q can be replaced by P 2 , the projection operator on £ 2 , 
and we get 

H' = — -^j- °<0|IH 1 P 2 IH 1 |0>° |O> 00 <O| + 0{k 3 ) . (2.18) 

It turns out that it is convenient to have a pictorial description of the processes in strong coupling 
perturbation theory. To evaluate the matrix elements one has to determine the action of the pertur- 
bation on the initial state on the r.h.s. and match the result with the final state on the l.h.s. In figs. 1 
the sucessive process of this evaluation is visualized by going along the vertical axis from the bottom 
to the top. Each vertical line represents a state where the spin is flipped down at a certain lattice 
site x G 7L d . If it is not affected by the perturbation the configuation remains unchanged. The corre- 
sponding diagram is depicted in fig. 1 (a). The hopping term ai(x)ai(x-\-k) flips the spins on adjacent 
sites, indicated by a horizontal line in fig. 1 (b) - (d). If it acts on an initial state where the spin is 
down at site x this spin is flipped up and his neighbour in the positive ^-direction is flipped down. 
Correspondingly the signal is moved along the negative ^-direction in case that the initial state has a 
spin down at site x+k. Finally, if both at x and x+k the spins are up initially they are flipped down by 
the action of ai(x)ai(x-\-k). Following these rules the graph contributing to (2.18) is easily found and 
depicted in fig. 2. The action of IH x on the initial vacuum creates an intermediate two-particle state, 
denoted by a dashed horizontal line. The spins are flipped back again by the subsequent action of the 
perturbation and we end up with the final vacuum state. Because there are dL d pairs of adjacent sites 
throughout the lattice, the graph carries a counting factor of dL d . We find 

H' = -K 2 jL d \tt}° °(n\ + 0{k 3 ) , (2.19) 

and the ground-state energy is 

E = -K 2 ^L d + 0(k 3 ) . (2.20) 
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Though E diverges as L — ► oo, this does not cause any trouble as long as the energy difference between 
the ground-state and the excited states stays finite. For later use we construct the next-to-leading 
order approximation of the ground state by means of Bloch's linear operator U (cf. eq. 2.16) 

d 

|0> = U |O>° = + -^2^2\ x ,x + k)° + 0{k 2 ) . (2.21) 

k = l x 

So when k is small but not zero the vacuum begins to mix with nearest-neighbour two-particle exci- 
tations. 

-k 



k 

Figure 2: Next-to-leading order contribution to the vacuum energy. 
2.4 The one-particle sector 

Next we consider the eigenspace Si spanned by the one-particle excitations | x G T} of energy 2. 

The projection operator is 

P^E^r^l , (2.22) 

X 

and the leading order hermitean reduced Hamiltonian reads 

H; (1) = J2 V|IHi|z 2 )V>° VI . (2.23) 

X 1 ,x 2 

The non- vanishing contributions are described by the elementary graphs figs. 1 (b) and (c). The initial 
flipped spin is moved one site in the positive or negative ^-direction, respectively, hence 

d 

H; (1) = (\x) 00 (x + k\ + \x + k) 00 (xfj . (2.24) 

k = l x 

In next order of perturbation theory the intermediate configurations are three-particle excitations of 
energy 6. Here the first contribution is made up of those graphs where the initial one-particle state 
la; 1 ) remains undisturbed while a vacuum fluctuation occurs on a distant pair of sites. This situation 
is shown in fig. 3 (a). Because there are d(L d — 2) distant sites the counting-factor is d(L d — 2). The 
second contributions come from those graphs where the flipped spin of the initial state is shifted over 
two lattice sites leading to the final states \x 2 ) = {x 1 ± k ± l)° and \x 2 ) = {x 1 ± k ^ They are 
represented by the graphs fig. 3 (b) - (e) with counting factor 1. We obtain 

h; (2) = -±(L d -A)?i-- A Yl E (\x} 00 (x+k+i\ + \x+k+i) 00 (x\ 

x k,l=l 

+ |a;) 00 (a; + A;-/| + |a; + A;-/) 00 (a;|) . (2.25) 

The additional contribution to the first term is due to the fact that for the graphs fig. 3 (d) and (e) 
the case k = I is excluded. For later convenience we prefer to write the summation over the directions 
without restriction and therefore have to add |c?Pi as a correction. 
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Figure 3: Next-to-leading order graphs with initial state {x 1 ) and final state \x 2 ) . 



The next step within Bloch's programme is the diagonalization of the reduced Hamiltonian 



H' = kH 



' (1) + k 2 H' (2) 



0(k 3 



(2.26) 



This is achieved straightforwardly if we pass to the coordinate space. Any state \ip)° £ Si has a unique 
expansion 

|^) = $>(z)|z>°, (2.27) 

X 

where ip(x) is a complex lattice function. The coordinate space representation of the Hamiltonian 
(2.26) is determined as follows 

(E[il>)(x) = "(zIH;^) . (2.28) 

To make the result more readable it is convenient to introduce forward and backward difference 
operators 

(VMx) = ^(x + k) - 1>{x) , (VW)(x) = i>(x) - i,(x-k) , (2.29) 
and the lattice Laplacian 



A 



def 



k ■ 



k = l 



In terms of these we find 



Hi = k{A - 2d} - K 2 \ ^A 2 - dA + d(d - 1) + ^L d \ + 0{k 3 ) . 



(2.30) 



(2.31) 



One recognizes that the vacuum energy (2.20) contributes to (2.31). Subtracting it from we define 
the resulting operator as a new effective one-particle Hamiltonian which no longer depends on the 
linear extent L of the lattice 

Hi = Hi - E . (2.32) 

Consequently, the infinite volume limit can be taken trivially and we can extend our analysis to the 
infinite lattice T = 7L d . In this case the static one-particle sector is represented by the Hilbert space 
of complex wave functions on 7L d satisfying 

\\n 2 = El^)! 2 <~ • (2-33) 

X 

Because Hi commutes with space translations by integer vectors (\J v ip)(x) = ip(x + v), its eigen- 
vectors are plane waves exp(ipa;). The momentum p ranges over the first Brillouin zone B = {p \ — it < 
Pk < }• Using the abbreviations 



P 



■2 def 



k = l 



E~ ^ ^ def . / Pk 

Pk Pk , Pk = 2 sin [ — 
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the eigenvalues read 

e^p) = k (p 2 - 2d) - ^ {(p 2 ) 2 - 4dp 2 + 4d(d-l)} + 0(k 3 ) . (2.34) 

The physical rest mass of the one-particle excitation is defined as the zero momentum energy 

m = £i(0) = 2 - 2dn - K 2 d(d - 1) + 0(k 3 ) . (2.35) 

At some critical inverse coupling k = k c (2.35) vanishes and the Ising model undergoes a continuous 
phase transition. In the following we will always assume that k < k c . In this case the one-particle 
energy is a positive, smooth periodic function of the momentum p. The kinetic mass to* of the particle 
is defined through the small momentum expansion 

2 

E^p) = m + + 0(p 4 ) . (2.36) 

2TO* 

In the present case one obtains 

to; 1 = 2k(1 + dn) + 0{n 3 ) . (2.37) 

At strong coupling the kinetic mass is much larger than the rest mass of the particle. This precisely 
reflects the fact that the one-particle configurations do not like to move, and in the limit k — ► 0, they 
are just static objects sitting on the sites of the lattice. 

To check the reliability of our method we discuss the physical properties of the resulting one-particle 
states. The translationally invariant static (improper) one-particle states read 

\p)° = (2ir)-iJ2 eipX \ x )° ■ ( 2 - 38 ) 

X 

They are normalized such that °{p\q)° = S(p — q) with S(p) denoting the ordinary (d- dimensional) 
Dirac ^-function. The next-to-leading order approximation of the true (improper) one-particle states 
is constructed as 

\ P ) = \p}° = \p}° + J-^p- J2 y, y+k)° + 0(k 2 ) . (2.39) 

k = l x,y 

Here the prime indicates the restriction x ^ y and x ^ y + k such that the summation is taken 
over three-particle configurations only. Below this notation is always assumed, i.e. in primed sums 
over static multi-particle states coinciding arguments are excluded. Normalizable one-particle states 
can then be formed by smearing with any smooth function f(p) which is periodic in all momentum 
components with period 2ir 

\f}= j dW(p)|p> = (2vr)^/( a; )| a; ) + K ^^^'/(^)k,y,y+A ; > + ^ 2 ) • (2.40) 

JB x k = l x,y 

Since the Fourier transform f(x) of f(p) falls off rapidly as \x\ — ► oo, one intuitively expects that the 
particle in this state is (essentially) localized in a bounded region 1Z of space which depends on the 
wave function f(p). The precise expression of this physical property is as follows (cf. ref. [16]). Let 
<D(x) be a local operator representing an 'observable' of the theory and assume that x is far outside the 
region 1Z. A measurement of this operator in the state |/) yields the value (/|0(a;)|/) in the average. 
Now if the particle is essentially confined to 1Z the measurement of <D(x) should give the same result 
as in the vacuum state Hence the one-particle state |/) £ 7i is well localized if 

(f\<D(x)\f) (f\f)(Sl\0(x)\Sl) (2.41) 

for \x\ — ► oo and any local observable <D(x). 

In appendix A we show that this condition is indeed satisfied. So we have good reason to argue 
that (2.40) represent the next-to-leading order approximation of the well localized one-particle states 
in the (d + l)-dimensional Ising model. 
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2.5 The two-particle sector 

For the determination of the two-particle scattering states we have to apply perturbation theory to 
the subspace of static two-particle excitations of energy 4 with projection operator 

P 2 = ^V^)°°(!/^I- (2-42) 

x,y 

To leading order the reduced two-particle Hamiltonian is given by 

H ' 2 (1) = iE' E' V^IIHil^TlsWH* 2 ,* 2 ! • (2-43) 

x 1 ,y 1 x 2 ,y 2 

The perturbation acts on the initial state la; 1 ,^ 1 ) by moving one of its flipped spins as described by 
the graphs in figs. 1 (b) and (c). The other excitation is not affected as represented by the elementary 
graph in fig. f (a). Their union are the graphs contributing to (2.43). Each of them carries a counting 
factor of 4 because a permutation x 1 <r+ y 1 in the initial state or x 2 <r+ y 2 in the final state leads to 
the same contribution. We find 

H' 2 (1) = "E'E {\x,y) 00 (y+k,x\ + \x,y+k} 00 (y,x\) . (2.44) 

x,y k = l 

The first contribution to the next-to-leading order correction is described by the graph in fig. 4 
(a). The initial state consists of two excitations on adjacent lattice sites y 1 = x 1 ± k, so the action 
of IH ! can flip both spins back. The final state is created from the vacuum at an arbitrarily shifted 
position. If the creation of the final state occurs before the initial state is annihilated this leads to an 
intermediate four-particle excitation depicted in fig. 4 (b). The permutation symmetry is accounted 
for by a counting factor of 4. Further contributions come from vacuum fluctuations (fig. 4 (c)) where 
the initial state remains unchanged. In case that its excitations are located on neigbouring sites the 
counting factor is 2(dL d — Ad + 1). Otherwise, the number of distant pairs of sites is decreased by 1 
and the counting factor is 2(dL d — 4d). Finally, the perturbation can shift one of the initial excitations 
over two lattice sites. The corresponding graphs are those of figs. 3 (b) - (e) completed by fig. 1 (a). 
Again each of them carries a counting factor of 4. We get 

H' 2 (2) = -^-8)P 2 -^^k,x + fc)°VMI + ^E E (\x,x + k) 00 (x + i,x\ 

x k = l x k,l=l 

~\~ | x , x k^j i^x ~\~ I ^ x | ~\~ | x , x ~\~ k^j i^x / , x | ~\~ | x , x k ^ ^ x / , x | ^ 

i / d 

-4E' E (\x,y) 00 (y + k + l,x\ + \x,y + k + i) 00 (y,x\ 
x,y k,l=l 

+ \x,y) 00 (y + k-i,x\ + \x,y + k-i) 00 (y,x\) . (2.45) 
Since P 2 is just the identity operator on £ 2 the infinite volume Hamiltonian is given by 

H' 2 = kH' 2 (1) + k 2 H' 2 (2) - E 2 P 2 + 0{k 3 ) . (2.46) 

As for the one-particle subspace, the diagonalization of the reduced Hamiltonian is most easily 
performed within the coordinate space representation. This means that £ 2 is identified with the space 
of normalizable wave functions on 7L d X 7L d which satisfy 

4>{x, y) = 4>{y, x ) , x ) = • ( 2 -47) 
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Figure 4: Next-to-leading order graphs with initial state la; 1 ,^ 1 ) and final state \x 2 ,y 2 ) . 

We extend the definition (2.29) of the difference operators by a subscript that indicates which argument 
of the wave function is shifted 

(^ X) k^)(. x i V) = ^(x + k, y) - ip(x, y) , (V^^X^ y) = ip(x, y+k) - ip(x, y) , (2.48) 
and correspondingly for V* k and V* k . Moreover it is convenient to introduce 

v„, M = v^v,,, + v„ iJfe + v,,, , v; M d ^ f -v; fc v;, + v; fc + v; ; , (2.49) 

where r] stands for x or y, respectively. Finally, for any r] £ 7L d we define the multiplication operator 



{Y)^){x,y) = 8(x-y-rj)il>(x,y) , 



(2.50) 



where 8(x) denotes the lattice ^-function which simply is Kronecker's 8 in case that the lattice spacing 
is equal to 1. 

Collecting everything together the reduced two-particle Hamiltonian reads 



H 2 = K (T 2 + W 2 ) = k {t^ + Wf* + k [T^ 2) + W^ 2) ] } + 0{k 3 



(2.51) 



where 



W 



(i) 



A, + A y - Ad , 
D A- fDoT^ + h.c.) , 



,(2) 



d A + A. 



2d(d-l) , 



W 



(2) 



rfD + (2d - 1) £(D £ + D_ £ ) - J D T^ - - [ D * ( V : 
k=i [ fc,;=i 



- V* 

y ,k ,1 x ,k ,1 



V 



x .k. — l 



V. 



y,k,-l 



} +h.c.| 



The first part of H 2 is just the sum of two effective one-particle Hamiltonians (2.32). Therefore T 2 
is identified with the kinectic term which determines the free propagation of two independent one- 
particle states. They interact through the potential W 2 which depends on the relative coordinate of 
the two-particles (x — y) and is rather short range. 

However, we note that due to the restriction for coinciding arguments, £ 2 ' 1S n °t the appropriate 
two-particle space. In particular, it is not possible to represent a state that describes two independent 
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one-particle wave functions moving freely. Instead, the natural Hilbert space for the evolution of 
a two-spin state is rather the symmetric tensor product of the one-particle subspace Si with itself. 
This space is represented by the symmetric wave functions without any further restriction and will be 
denoted by 7i 2 . Since S 2 is a closed subspace, 7i 2 is the direct sum of S 2 and its orthogonal complement 
S 2 , namely the space of symmetric wave functions which vanish except for coinciding arguments. 

Now the discussion of scattering is much easier if the free and the interacting dynamics are defined 
on the same Hilbert space. Otherwise one would have to use the formalism of two Hilbert space 
scattering theory as described in ref. [18], for example. Though conceptually satisfactory this method 
bears some complications in practical calculations. It is therefore convenient to extend H' 2 to a self- 
adjoint operator H 2 on 7i 2 . Since D projects on S 2 such an operator is given by 

H 2 = (]I-D )H' 2 (]I-Do) + kADo , (2.52) 

where A is an arbitrary real constant. Since the free Hamiltonian T 2 trivially extends to an operator 
on TC 2 we can rewrite (2.52) as 

H 2 = k(T 2 + V 2 ), (2.53) 
with a modified but still very short range potential 

V 2 = -D T 2 - T 2 D + D T 2 D + (11 - D )W 2 (1I - D ) + AD . (2.54) 

The wave functions ip £ S 2 are rapidly decaying. Consequently, they do not belong to the absolutely 
continuous spectrum but represent bound states with energy kX. Therefore all physical scattering 
quantities must eventually be independent of A. 

Here we make a technical remark. Since the restriction to the symmetric wave functions complicates 
the diagonalization of the Hamiltonian (2.53), we will henceforth consider H 2 as an operator on the 
entire space of normalizable wave functions L 2 (7L d X 7L d ). This is possible because H 2 is identical 
on 7i 2 and L 2 (7L d X 7Z d ), and commutes with the unitary operator Ui 2 representing the permutation 
(V i 2 ip)(xi, x 2 ) = ip(x 2 ,Xi). So after we have succeeded in finding the eigenstates of the Hamiltonian 
(2.53) defined on L 2 (7L d X 7L d ) the true two-particle sector of the Ising model is restored by projecting 
on the subspace of symmetric wave functions. 

Since the two-body Hamiltonian (2.51) commutes with the translation operator (\j v ip)(x,y) = 
ip(x + v,y + v) the center of mass motion can be separated. Let p x and p y denote the momenta 
corresponding to x and y, respectively. Introducing 

r = x-y, q=^(p x -p y ) and R = x + y , P = ^(p x + p y ) , (2.55) 

the simultaneous (improper) eigenstates of H 2 and U„ have the general form 

^(x,y) = e lPR ^(r) , (2.56) 

where ip £ L 2 (7L d ). Now if we introduce 

H(P) = e lPi? H 2 e- lPi? , (2.57) 

the action of H 2 is expressed as 

(K 2 f)(x, y) = e lPR (H(P)f) (r) . (2.58) 

So the eigenvalue problem of the two-particle Hamiltonian H 2 is reduced to the diagonalization of 
the effective one-body Hamiltonian H(P). The corresponding Hilbert space 7i(P) is a translationally 
invariant subspace of L 2 (7L d X 7L d ). It is identified with the space of normalizable wave functions in 
the relative coordinate r. For simplicity we will restrict our attention to the center of mass system, 
i.e. in the following we assume that the total momentum P of the two particles is zero. 
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The explicit calculation of (2.57) is straightforward. We state the final result 

H = k(T + V) (2.59) 

where the free Hamiltonian reads 

T = 2{A r - 2d} - | {A r 2 - 4dA r + 4d(d - 1)} + 0(k 2 ) , (2.60) 

and the potential is given by 

V = V (1) + kV (2) + 0(k 2 ) , (2.61) 



= -D T-TDo + D TDo + ADo , (2.62) 
V (2) = (2d-l)£{D s + D_ s } 

k 

1 d 

"2 T, { D * V *,M " D -* V r ,M + D ^V;,,,_ ; - D_jV r|ij _, + h.c.} . (2.63) 



3 The scattering formalism 

The two-particle dynamics of the (d + l)-dimensional Ising model is determined by the effective 
Hamiltonian (2.59). To obtain the desired strong coupling expansion for the elastic scattering states 
and the transition matrix we proceed following the rules of ordinary quantum mechanics. 

3.1 Brief review of quantum mechanical scattering theory 

The starting point of scattering theory is the physical assumption that the constituents in the initial 
and final state are so far away from each other that the interaction between them is negligible. In our 
case this means that the scattering states of the effective one-body system (2.59) are those vectors in 
the Hilbert space which are 'asymptotically' free in the distant past and future. The mathematically 
precise meaning of this statement is that for any / in the absolutely continuous subspace TC ac (T) of 
the free Hamiltonian T there is a vector f m £ L 2 (7L d ) such that 

li m II e" lT 7 - e" lH 7 m 11 = 0, 

t^ — co 

and similiarly for pairs of states {/, / out } that approach each other as / — ► +oo. Let P ac (T) denote 
the projection operator on TC ac (T). Then the basic existence question of scattering is reduced to the 
problem of proving that the M0ller operators 

fi in (H,T) = s-lime lHt e- m P ac (T) , (3.1) 

t— ► — CO 

fiout(H,T) = s-lime lHt e- m P ac (T) (3.2) 
exist. In case they do, in- and out-going scattering states are denned as 

r = n m f, f out = n out f. (3.3) 

Furthermore we expect that the ranges "H; n /out of the in- and out-going M0ller operators are equal to 
the absolutely continuous subspace 7^ ac (H) because any interacting state which looked like a free state 
in the distant past should again look like a free state for large positive times. This property is called 
the completeness of the M0ller operators. 
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To prove that 0; n and O out exist and are complete we take advantage of the special structure 
of the potential. From the derivation of the reduced two-particle Hamiltonian by means of Bloch's 
perturbation theory it should be clear that to all orders of k the effective potential V will be a sum of 
operators each of which is the product of a constant or a difference operator and a local multiplication 
operator D^. Now any difference operator on L 2 (7L d ) is bounded. Furthermore, is trace class 
because it is a positive operator and satisfies 

oo oo oo 

v — 1 v — 1 r v — 1 

for any orthonormal basis {(f v }™-i of L 2 (7L d ). The product of a bounded operator and a trace class 
operator is trace class and the family of trace class operators is a vector space (cf. theorem VI. 19 of 
ref. [18]). It follows that the potential V is trace class, and the prerequisites of the Kato-Rosenblum 
theorem are fulfilled. This guarantees the existence and completeness of fi; n / out (H, T) (cf. theorem 
XI.8 of ref. [18]). 

The fundamental operator of scattering is 

S = f ^out (^m) t : Wi n — ► Wout • (3.4) 

It maps in-going states onto out-going states with the same probability distribution of momentum. 
Hence the scattering operator correlates the past and the future asymptotics and determines the 
physics of the scattering process. The scattering operator commutes with the Hamiltonian H and is 
unitary on 7^ ac (H). 

The free Hamiltonian T has the purely continuous spectrum 

Ke 2 (q) = ti{Px) + ti(Py) (3.5) 

where 

e 2 (q) = 2{q 2 - 2d} - | {(g 2 ) 2 - Adf + Ad(d - 1)} + 0(k 2 ) . (3.6) 

From [15] we know that the series representation for T is convergent if the inverse coupling is sufficently 
small. It follows that the two-particle energy is a smooth function of the relative momentum q. The 
plane waves 

HT,r)=-^-r , -7T < q< 7T (3.7) 
{2ir) 2 

are a complete set of (improper) eigenfunctions normalized such that (<J>(g'), <J>(g)) = S(q' — q). In the 
sense of distributions (improper) in- and out-going scattering states can be defined as 

$ in (g) d ^ f ft m <%) , <P ut (g) d ^ f ft out <%) . (3.8) 

They are (improper) eigenstates of the Hamiltonian H with energy Ke 2 (q). They span the absolutely 
continuous subspace 7^ ac (H) and satisfy 

($ in (?'),* in (?)) = ($ ou V),* out (?)) = (<%'),<%)) • 

Starting from (3.8) one deduces an integral equation for the scattering solutions, the famous Lippmann- 
Sch winger equations 

$ m/out/ n = $ / n _ Um L W n/out (g) . (3.9) 

p\o T-£ 2 (?)Tip 

With the solutions of these equations at hand normalizable scattering states are obtained by smearing 
with a smooth periodic function in the relative momentum q 

f m = I d d qf(q)^(q) , / out = / d d qf(q)^(q) • 
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This leads to an integral representation of the scattering operator 

(r,sr) = / ay / d d q f( q rf(q)^ m (q'),s^( q )) , 

JB Jb 

and the physics of the scattering process is determined by its kernel, the S-matrix 

S(q',q) = ($ m (y),S$ m (g)) = ($ m (g'),$ out (g)) . (3.10) 

By means of the Lippman-Sch winger equations one can establish an explicit formula 

S(q<, q) = ($(?'), <%)) " 2vri« S[E 2 (q') - E 2 (q)]T(q', q) , (3.11) 

where E 2 (q) is the full two-particle energy 

E 2 (q) = 4 + Ke 2 (q) , (3.12) 

and the on-shell transition matrix is defined as 

T(q',q)= (<%'), V<I> ln (g)) . (3.13) 

Finally we have to remember that only the symmetric wave functions belong to the true two- 
particle sector of our model. Since the permutation operator Ui 2 commutes with both T and H, it 
also commutes with the M0ller operators (3.1), (3.2) and hence with the scattering operator (3.4). This 
means that in (3.11) and (3.13) we simply have to substitute the wave functions by their symmetric 
counterparts 

x def 1 r / x / s-, COs(or) 

**(q;r) d = f -{<%;r)+<%;-r)} = — ^ , (3.14) 
K(q,r) = ^{^(q;r)+^(q;-r)} , (3.15) 



which satisfy 



:$ 8 ( g '), *s(?)) = (<W), K(q)) = \{Hq'-q) + %'+?)} • (3.16) 



3.2 Generalized Born series 

It is convenient to introduce the Green's operator 

G(q) = lim . . 1 m . (3.17) 

KHJ p\o e 2 (q)-T + ip V ; 

Since the free Hamiltonian T commutes with the unitary operator ~Ur representing the hypercubic 
transformations 

(tWO(r) = ^(Rr) , R G 0(d, TL) , (3.18) 

so does G(q). From now on we follow the usual convention that p > is understood as the limit 
p 0. In coordinate space (3.17) is referred to as the Green's function 

G(q;r,r') = {r\G(q)\rr (3.19) 

For fixed q it is a well defined function on (7L d X 7L d ) which is invariant under identical hypercubic 
transformations in both arguments and has an integral representation as 

r (\ d r> pip(r-r') 

G(q;r,r')= / 7^—-^ (3.20) 

Vy ' ; Jb (2vr) d e 2 (q)-e 2 (p) + ip V ; 
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Since the energy e 2 (q) is a smooth hypersurface in momentum space the Green's function is a smooth 
function of the relative momentum q, at least as long as the energy hypersurface does not intersect with 
the boundary of the Brillouin zone. For notational convenience we write G(q; r) instead of G(q; r, 0). 
Using (3.17) the Lippmann-Sch winger equation for the in-going scattering states reads 

$ in ( g ) = $( g ) + G( 9 )V$ in ( 9 ) . (3.21) 

The potential is a convergent power series in the inverse coupling where the range of V increases with 
increasing order in k. We are heading now for a solution of the above integral equation 1 . The first 
step is to decompose (3.21) as follows 

$ in ( g ) = + G(q) V (1) $ in ( 9 ) + G(q) [V - V (1) ] $ in ( 9 ) . (3.22) 

Since for the leading order potential = 0(1) the second term is of the same order as the free 
solution $(q), whereas the third contribution is proportional to k. In the next subsection we show 
that 

f;"( g ) = f(g) + G(g)V«f;"( g ) (3.23) 
can be solved exactly thanks to the locality of V^K Substituting this result, eq. (3.22) becomes 

$ in ( g ) = $f( q ) + []| - G(g)VW] _1 G(q) [V - V^] $ in ( 9 ) . (3.24) 

Thereby we have derived a modified integral equation for the scattering solution. To solve eq. (3.24) 
we make a Born series ansatz 

$ in (?) = £ [l- G(g)VW] " G(q) [v - V^] t>?(q) (3.25) 

which represents a power series in the inverse coupling. For a discussion of convergence we consider 
(3.25) in the coordinate space representation. Since all contributions to the potential V have finite 
support there are no infinite sums but products of Green's functions evaluated at definite coordinate 
points. To any order 0(K n ) these products consists of n factors at most. In addition the Green's 
function is bounded. So there is no doubt that the Born series converges if the inverse coupling is 
sufficiently small. Consequently, (3.25) is indeed a solution of the Lippmann-Sch winger equation (3.21) 
and a valid representation of the Ising model's scattering states. 

3.3 Leading order 

We consider the leading order Lippmann-Sch winger equation (3.23) in the coordinate space represen- 
tation 

^T(q; r) = <%; r) + £ G(q; r, r')V (1 ^T(q; r') . (3.26) 

r 

To evaluate the r.h.s. we substitute the explicit form of the leading order potential (2.62). By means 
of the integral representation (3.20) one verifies the difference equation 

{-T + e2 (q)}G(q;r) = S(r) . (3.27) 

for the Green's function. It follows that 

J2 G(q; r, r<) T D *™(q; r<) = z[ x) {e 2 (q)G(q; r) - 6(r)} . (3.28) 

r' 

1 From what follows it will be clear that we can derive a solution for the out-going states analogously. 
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where we have abbreviated = ^"(g; 0) Furthermore, inserting a complete set of plane waves we 
obtain 

d d p 



Do T Do $?"(?; r) = Ij^rMP) *P %) 



3.29) 



Finally the action of D T on the wave function $™(q;r) leads to a a sum of values of $™(q;r) at 
definite coordinate points r which we abbreviate by a complex number z~ 2 ^ 

D T^ n (g;r) = 4 1 ^(r) . (3.30) 

Substituting (3.28)-(3.30), eq. (3.26) becomes 

r) = <%; r) + z^ 6(r) - {z^ + [e 2 (q) - u - A] z^} G(q; r) . (3.31) 

So the coordinate dependence of the leading order scattering states is determined by the free solution 
<i>(g; r) and the Green's function G(q;r). To determine the values of the unknown numbers and 
z~ 2 ^ we first evaluate (3.31) at r = 

{e 2 (q) -u- A}G(g; 0) z ( p + G{q; 0) z { 2 V) = 0) , (3.32) 

and then act on (3.31) by D T 

{A - e 2 (q) + [e 2 (q) -u- \]e 2 (q)G(q; 0)} z ( p + e 2 (q)G(q- 0) z { 2 V) = e 2 (q)<S>(q- 0) . (3.33) 

Combining (3.32) and (3.33) one obtains a simple matrix equation for and z~ 2 ^ which can be solved 
easily 

4 1} = 0, 4 1} = b^GM)!" 1 . (3.34) 



In particular we find that the scattering solutions vanishes at r = 0. We substitute (3.34) into (3.31) 
and symmetrize the result. The leading order scattering solution reads 

r) = $ s (g; r) - G ^ . (3.35) 

With this solution at hand it is also straightforward to calculate the leading order on-shell transition 
matrix (cf. eq (3.13) 

T 1 (q',q)= ($,(?'), V^S^g)) = - [(270^(9; 0)]" 1 . (3.36) 

The reader should not be frightened at the apparent violation of time reflection symmetry q 1 ^ q. 
From (3.20) it follows that G(q'; r) = G{q, r) if e 2 (q r ) = e 2 (q). 

3.4 Next-to-leading order 

The next-to-leading order contribution to the Born series (3.25) reads 

= [a - G(g)VW] _1 G(q) [V - VW] $T(q) . (3.37) 
In coordinate space we obtain the integral equation 

r) = nA(q- r) + J2 G(q; r, r')V^ 2 n (q; r<) (3.38) 

where we have abbreviated 

A(q- r) = Y, G(q; r, r')V (2 ^T(q; r') , (3.39) 

r' 
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and the next-to-leading order potential is given by (2.63). Except for the inhomogeneity, eq. (3.38) 
has exactly the same structure as the leading order Lippmann-Sch winger equation (3.26), hence 

r) = KA(q; r) + z^ 6(r) - {z™ + [e 2 (q) - u - X] z[ 2) } G(q; r) , (3.40) 

where 

z[ 2) = ^(q;0), D T^(?;r) = 4 2) %). (3.41) 
Substituting the leading order solution (3.35) and the potential (2.63) it is straightforward to calculate 



r ■ 



■) d/2 r G(q,0) } 

-(^7?E sin (ft){G'(g;r-fc)-G'(g;r + A ; )} . (3.42) 

Here and below we make use of the identity G(q; ±A;) = G(q; 1) which is a consequence of the hypercubic 
symmetry of the Green's function. 

Since the coordinate dependence of the inhomogeneity A(q; r) is again determined by the Green's 
function, eq. (3.40) can be solved exactly in the same way as (3.31). Using (3.27) one finds 

D T A(q; r) = e 2 (q)A(q; 0) 6(r) . (3.43) 

So the resulting equations that fix the unknown variables z^ and z^ are the same as (3.32) and 
(3.33) if we only replace <I>(g; 0) by KA(q; 0). It follows 

(2) n ( 2 ) A(q] °) fo AA\ 

A =0 ' 4 =k gm)' (3 - 44) 

Substituting these values into (3.40) and symmetrizing the result we end up with 

r) = « [Mi; r) - |=^°2 G(q; r)} , (3.45) 



where 



^ (9 ' r) = -^{^" 2d+1? ^^}?{ G(9;r -* ) + G(9;r + 

^— £cos( 9ife ) {(?(?;/•- £;) + G(g; r + £;)} , (3.46) 



(2vr) d / 

(2d-i)G( g; i) OJ , OJ g(g;i) l 

The next-to-leading order contribution to the transition matrix consists of two parts 

T 2 (g', q) = ($,(?'), V^*^?)) + « ($ s (g'), V (2) $^(g)) . (3.48) 
Since <j> 2 n s (g; 0) = the first term reduces to 

TfV, ?) = - (Ml'), D T$^ 2 (g)) . (3.49) 
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which is calculated trivially using the definition of z 2 (3.41) and the result (3.44). To determine 
the second contribution we simply have to insert the result for V^^^g; r) obtained previously 
(cf. eqs. (3.39) and (3.42)). Combining everything together the final result reads 



f2vr) 



2^cos(g fc )cos(^J 



(y a -2 < o(S B -2 < o- (M - 1)G(9;1) 



G(q;0) 



-2.-2 , , . ojGt?; 1) 

q + q — Ad + 2d- 



G(q;0) 



(3.50) 



The approximation for the scattering quantities presented here relies on the Born series repre- 
sentation of the scattering solution (3.25). It is strongly motivated by the physical structure of the 
effective one-body quantum system derived by means of Bloch's perturbation theory. Alternatively 
one could think of a polynomial approximation in the inverse coupling which would require to expand 
the Green's function (3.20) in powers of k. This is perfectly alright here but demands some care if 
G(q; r) is more complicated. It might even happen that the Green's function is singular at k = 0. 
Anyhow, a polynomial approximation for the scattering quantities is less natural and will be worse in 
general. 



4 The (1 + l)-dimensional model 

In the preceding section we have derived the general formulae for the leading and next-to-leading 
order contributions to the scattering solution and the transition matrix. To get a feeling of the results 
obtained we specialize to the simplest case of d = 1 spatial dimension. For this model exact solutions 
exist [19] and it is therefore a good testing ground for perturbative methods. 

The next-to-leading order approximation for the two-particle energy (3.6) can be written as 

E 2 (q) = 4 - 4k cos g + 2k 2 sin 2 q . (4.1) 

As long as the inverse coupling is smaller than the critical value k c = 1 the denominator of the Green's 
function (3.20) factorizes 

1 r dv e ipr 
G(q;r)=— —- — . (4.2) 

2k J-k 2tt (cosp — cos q + ip)[cosp + cos q + 2/k) 

For fixed q and r the integrand extends to an analytic function of p in the domain C \ {pi,p 2 } where 

p 1 = \q\ . ^ ^ (mod 2ir) , p 2 = vr + i arcosh(cos q + 2/k) (mod 2ir) . 

By means of the fundamental theorem in the calculus of residues one obtains 



-l e 



ikIM 



(-l) r e-' , M 



G(q;r)=— < — — — — ik- '- \. (4.3) 

4(l + Kcosg) I sin \q\ fj, \ 

where we have abbreviated 

H = arcosh(cos q + 2/k) , p, = k sinh(/i) = ^4(1 + k cos q) — k 2 sin 2 q . (4-4) 

With this result at hand it is now an easy task to determine the next-to-leading order approximation 
of the scattering solution. Using (3.35) and (3.45) we find 

(q; r) = |cos(gr) + Cl e l|3l|r| + c 2 (-l) r e -' ,|r| ) + C(k 2 ) , (4.5) 



/2vr 
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where the coefficients C\ and c 2 are given by 

-f 



Cl 



c 2 



H — ik sm \q\ 

\k sin \q\ 
p, — \k sin \q\ 



2\k sin \q\(l + k cos g) 
jl — ik sin |g| 

2ik sin |g|(l + k cos q) 
jl — ik sin |g| 



(4.6) 



In order to calculate the phase shifts we have to diagonalize the S-matrix (3. II). Since the scat- 
tering operator commutes with the two-particle Hamiltonian it is favourable to introduce energy 
eigenstates \E 2 ) and \E 2 in) which are normalized such that 



{E' 2 | E 2 ) = (E' 2 in | E 2 in) = S(E' 2 - E 2 



(4.7) 



The two-particle energy (4.1) is a even one- dimensional function of the relative momentum q. Conse- 
quently it is E 2 = E 2 if and only if q 1 = ±g, hence 



HK ~ E 2 



dE 2 



dq 



{S(q< -q) + S(q< + q)} . 



(4.8) 



It follows that the desired eigenstates are obtained if we define 



\E, 



ldE 2 



2 dq 



-1/2 



$ s (g) , \E 2 in) 



ldE 2 



2 dq 



-1/2 



(J) 11 



(4.9) 



With respect to these new bases of free and scattering states the scattering matrix takes the form 



S(E' 2 , E 2 



1 — 47HK 



dEo 



dq 



T(E 2 ) \ S(E' 2 - E 2 



(4.10) 



Using the general results (3.36) and (3.50) for the transition matrix we obtain the next-to-leading 
order approximation 



S(EL E 2 



4(1 + k cos q)(l — \k sin |g|) 
(jl — ik sin |g|) 2 



S(E' 2 - E 2 ) + 0(k 2 ) 



(4.11) 



Due to the unitarity of the scattering matrix its eigenvalue has absolute value 1. Therefore it can be 
parametrized as exp(2i<*)) where the angle S is referred to as the phase shift or the scattering phase, 
respectively. For any relative momentum q and any inverse coupling k < k c eq. (4.11) allows an 
approximate calculation of the phase shift in the (1 + l)-dimensional Ising model. Within errors we 
find the constant value 

S = - • (4.12) 



7T 

2 



In this simple model it is worth while to consider the polynomial approximation as well. The result 



is 



[q;r) 



/2vr 



for the scattering solution and 



6 



7T 

2 



sm(\q\\r\) + 0(k 2 ) 



0(k 2 



(4.13) 



(4.14) 



for the phase shift. In fact (4.13) and (4.14) hold not only in the center of mass system but for any 
total momentum P of the scattering particles. An interesting observation is of course that there is 
no O(k) contribution. This result may be compared to the strong coupling series for the one-particle 
mass. From our general expression (2.35) it follows that for d = 1 there is no second order contribution 
either. Moreover, higher order calculations show that the series actually truncates after the first term 
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and one finds that the leading order contribution represents the exact result for the one-particle mass 
[21] 

m = 2(1 - k) . (4.15) 

It may be that this is also the case for the scattering matrix meaning that —1 is indeed the exact 
result in the symmetric phase k < k c = 1. The assumption is corroborated by recent studies of Lang 
et al. [7]. These authors determined the phase shift to be equal to 7r/2 by means of Liischer's finite 
size technique approach [1] alluded to in the introduction. Anyhow, in ref. [22] it was shown that at 
the phase transition point k c = 1 the scattering matrix is —1. 

We end our investigations of the (1 + l)-dimensional model with a discussion of the physical 
properties of the derived scattering states. The (improper) eigenstates of the hermitian reduced two- 
particle Hamiltonian read 



\P*,P 9 m)° = J2 elP(X+y) sin (MI* " V\)\x, y}° + 0(k 2 ) . (4.16) 

They are labelled by the momenta p x and p y of the colliding particles and satisfy 

°(p' x ,p' y m\ Ps , Pv in) = °(p' x \p x } 00 (p' y \p y }° + °(p' x \p y } 00 (p' y \p x }° ■ 

The corresponding basis of true in-going scattering states is constructed by means of Bloch's linear 
operator U 2 and found to be 

\p x ,p y in) = \p x ,p y in) + ^ i^S(p x + p y )\tt)° 

x,y,z J 

In terms of these, normalizable in-going scattering states can be written as 

&Px I dp y f 1 (p x )f 2 (p y )\p x ,Py in) . (4.17) 

- 7T J — 7T 

where the one-particle wave functions fi and f 2 are supported in momentum space in such a way that 
the corresponding ranges of group velocities dEi(p)/dp have opposite sign. The physical characteri- 
zation of |/i,/2 in) is that for large negative times it evloves into a state which describes two widely 
seperated one-particle wave packets approaching each other. By an argumentation similar to the one 
applied in the case of localized one-particle states we expect that for any local observable ®(x) the 
in-going scattering states \ fi, f 2 in) should satisfy (cf. ref. [16]) 

(A,/ 2 in|0(x)|A,/ 2 in) = (AlAX/slOC^)!^) + {/ 2 |/ 2 >(/i|©(^)|A> 

-(fi\fi)(f2\f2)(n\0(x)\n) , (4.18) 

up to terms which vanish rapidly when the distance between the wave packets fi and f 2 is made large. 
In appendix C it is shown that the states (4.17) indeed fullfil this basic physical property. 



5 The (3 + l)-dimensional model 

The difficulty with the higher dimensional models is that we do not have an analytical solution 
for the Green's function (3.20). Therefore we will not be able to present the next-to-leading order 
approximation of the in-going scattering states in a closed form as for the simple (1 + l)-dimensional 
case. However, due to the locality of the potential, to calculate the leading and next-to-leading order 
transition matrix the Green's function is needed only at a few points near the origin. 
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5.1 Small momentum expansion 

Bloch's perturbation theory yields the two-particle energy E 2 (q) as a power series in k which is 
convergent if the inverse coupling is sufficiently small. We introduce momentum spherical variables 

q 1 = q sin(0) cos(^) , q 2 = q sin(0) sin(^) , q 3 = qcos(0) . (5-f) 

Since the momentum cartesian coordinates are restricted to the interval q k £ ( — 7r, x] , for q > it the 
angular variables do not cover the whole range < 9 < it and < (f> < 2ir, respectively. Instead 
they must be restricted to some subdomain of the surface of the sphere since, pictorially speaking, 
the sphere of radius q > it intersects the cube of edgelength 2ir. Physically the range q > it is the 
region where the model is dominated by lattice artefacts. But if the abolute value of the relative 
momentum is less than 7r the two-particle energy is a smooth function of q. It follows that E 2 (q) 
has a convergent series representation in powers of q. In the limit of small momentum the continuum 
rotational symmetry is restored on the lattice. Therefore the first non-constant contributions to this 
series is proportional to q 2 (cf. eq. (2.36)), hence 

2 00 

E 2 (q) = 2m+ + Y J h(n;K)q» . (5.2) 

m * u=3 

It is convenient to intoduce a modified absolute value of momentum 

q* = f \fm* {E 2 (q) - 2m} . (5.3) 

At low energy or small momentum, respectively, one has q* = q, so (5.3) is a suitable parameter for 
on-shell small momentum expansions of the scattering quantities. From the implicit function theorem 
we conclude that eq. 5.2) can be inverted for small values of q*. This leads to a series expansion for 
the (true) absolute value of the relative momentum of the scattering particles 

oo 

1 = 1* + c v (fi|K)?* • (5.4) 

v = 2 

The coefficients c„ are polynomials of sin(#), sin(^) and cos(#), cos(^). Therefore they can be expressed 
as sums of spherical harmonics Y, m (fi). 

In order to diagonalize the scattering matrix one proceeds as in the (1 + l)-dimensional case. We 
define free states \E 2 ,£l) and scattering states \E 2 ,£l in) through 

\E 2 ,n) = q(j^ V2 $„(«,), \E 2 ,nin) = q(j^ ^ *?(q) . (5.5) 



From the normalization (3.16) it follows that 



-1/2 



(E' 2 ,n'\E 2 ,9,<j>) = (E 2 ,n' in \E 2 ,9,<j> in) = (^^J %' - g) tf s (fi' - fi) , (5.6) 
where 

W - fi) = ^r^y {6(9'-9)6(4>'-4>) + 6(0' + - ir)6(<f>'-<t> - vr)} . (5.7) 

The r.h.s. of eq. (5.7) is zero unless the primed angular variables either coincide with the unprimed 

ones or are related to them by a space reflection. In both cases E 2 = E 2 if and only if q 1 = q. So can 
substitute 

%' ~q) = ~ E 2 ) (5.8) 
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on the r.h.s. of eq. (5.6) and end up with 



(E' 2 , ft' | E 2 , ft) = S(E' 2 - E 2 ) S s (tt' - ft) . (5.9) 
With respect to the new basis of states the scattering matrix is given by 

S(£ 2 , ft'; E 2 , ft) = S(E 2l ft', tt)S(E' 2 - E 2 ) , (5.10) 

where 

/A TP A TP \ — 1/2 

S(£ 2 ;ft',ft) = 6 a (Sl' - ft) - 2irmqq' {jj^-^ J T(£ 2 ; ft', ft) . (5.11) 

Remembering the results of section 3 we known that the transition matrix is determined by the 
Green's function at definite coordinate points plus some kinematical contributions. Furthermore, the 
integral representation (3.20) together with (5.2) allows us to expand G(q;r) as a convergent power 
series in q if the relative momentum is sufficiently small. By means of (5.4) one then deduces a series 
representation for the scattering matrix at fixed energy 

00 

S(E 2 ; ft', ft) = 6 s (tt - ft) + J2 0) € ■ (5-12) 

The coefficients S v are sums of spherical harmonics in the momentum angular variables (ft', ft), where 
the maximum angular momentum quantum number / increases with increasing order in q* . Eq. (5.12) is 
a general result independent of the order of truncation for the effective Hamiltonians. The calculation 
of the first few orders provides an insight into the small momentum behaviour of the model's scattering 
quantities. For instance, one will be able to determine the scattering length and get an idea of the 
angular dependence of the scattering matrix. 

To be explicit we perform the small momentum expansion up to order O(ql) using the results from 
the next-to-leading order approximations derived in the preceeding sections. The two-particle energy 
reads 

E 2 (q) = 2m+^-^-(q 2 f , (5.13) 

to* 2 

where 

to = 2-6k(1 + k), to; 1 = 2k(1 + 3k) . (5.14) 

It is a specialty of this approximation that E 2 is a polynomial in q = \J~if ■ In appendix C we argue 
that the Green's function has a convergent series representation 

00 

G(q-r) = Y,G v (r)q (5.15) 

and calculate the coefficients up to the fifth order to 10 significant digits. The results for r = 0, r = 1 
and some values of k are listed in tables 2 and 3 of appendix C. Using the definition (5.3) we obtain 
from (5.13) 

q = q, + -K 2 m* ql + ql + 0(ql) . (5.16) 

Together with (5.16) eq. (5.15) provides a small momentum expansion for the Green's function. To 
handle the remaining terms in eq. (5.11) we need to expand the absolute value of the relative momen- 
tum (cf. eq. (5.4)) 

q = q* - ^b 4 m* ql + (^-b\ml - ^b 6 m^j ql + 0(ql) , (5-17) 

where 

h = -y - f (1 + 3K)/i(fi) , h = ^^(ft) + ^(1 + 3*)/ 2 (n) , 
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and 

= sin 4 cos 4 <j> + sin 4 sin 4 <j> + cos 4 , (5.18) 
/ 2 (fi) = sin 6 cos 6 + sin 6 sin 6 + cos 6 . (5.19) 

Combining the above results the desired low energy expansion for the scattering matrix is now 
obtained by some lengthy but straightforward algebra. We state the final result 

5 

S(E 2 ; fi', fi) = 6 s (tt -tt) + J2S„ ql + O(qt) , (5.20) 



(5.21) 



and 



where 

f _ «11 c _ «31 + «32[/l(» / ) + /l(»)] 

"Jl — 9 7 ^3 — 9 7 

c _ «21 c _ «41 + «42[/l(^) + /l(»)] 

^2 — 5" 7 ^4 — ^ 7 

5*5 = {« 5 i I a 52 [/i(fi') + /i(fi)] + asal/^') + A 2 (ft)] + a M /i(ft')/i(fi) 

+a 55 [/ 2 (fi') + / 2 (fi)] + a 56 / 3 (fi', fi)} • (5.22) 

The function / 3 appears when expanding the angular dependent part of the next-to-leading order 
transition matrix (3.50) 

/ 3 (fi', ft) = sin 2 0' sin 2 cos 2 <j> cos 2 0' + sin 2 sin 2 0' sin 2 <j> sin 2 0' + cos 2 cos 2 0' . (5.23) 

The coefficients a^- are functions of k and Gr„(r), for instance 

. m« Co(0)-30k{G'o(0)-G' 1 (0)} 2 

ail = 1 T Gjoy • (5 - 24) 

The higher order contributions become more and more lengthy and are unsuitable for publication. If 
one substitues a certain value for the inverse coupling together with the results of tables 2 and 3 of 
appendix C they are just complex numbers. 

For the diagonalization of the scattering matrix we expand (5.20) into spherical harmonics. Since 
we are confined to the subspace of symmetric wave functions we have to restrict ourselves to even / 
values. Hence 

oo 21' oo 21 

S(£ 2 ;fi',fi) = £ J2 E E Ytf(SiyS(E 2 ) 2l , m , t2lm Y£(Sl) , (5.25) 

('=0 m' = -2l' 1=0 m = 2l 



with 



and 



S(E 2 ) 2V m '_ 21m — du'fimm' + E m',2f m 1* + C(?*) 7 (5.26) 



«^ m , 2/m = / dfi' / dfty™ (n')s v Y™(ny . (5.27) 



Up to the second order (5.20) does not depend on the the angular variables. So the correction to 
the identity matrix in (5.26) is only for the first eigenvalue. As a consequence the scattering length 
can be determined straightforwardly 

ao = i im ^W>° = ^ii . (5.28) 
3.-0 2i ivr y J 
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2 4 6 8 10 

k X 100 

Figure 5: The geometric cross section as a function of the inverse coupling. 

Fig. 1 shows the geometric cross section a = Aira\ as a function of the inverse coupling k. From (2.35) 
one obtains k c = 0.264 + 0(k 3 ) for the (3 + l)-dimensional model. Therefore it is not reasonable to 
go beyond k = 0.1 in a strong coupling expansion. The interaction strength increases almost linearly 
with increasing magnitude of k. Nevertheless the sensitivity is not very strong since an increase of a 
factor of 10 only results in a change of 11 per cent for the cross section. The corresponding area is 
about 1.3 plaquettes on the lattice. 

To evaluate the higher order contributions (5.27) we write the functions as sums of spherical 
harmonics and use the orthogonality relations. To third and fourth order the maximum value is / = 4 
leading to a non-trivial (15 X 15) submatrix. If we include the fifth order contribution the maximum 
value increases to / = 8 and we obtain a (45 X 45) submatrix. 

5.2 The phase shifts 

The scattering matrix is invariant with respect to cubic transformations. As a consequence its eigen- 
states come as multiplets which transform according to the irreducible representations of the cubic 
group 0(3, 7L). It is well known that the special cubic group SO(3, 7L) has 5 irreducible representations. 
They are denoted by Ai,A 2 ,E,Ti and T 2 , and their dimensionalities are 1,1,2,3 and 3, respectively. 
The irreducible representations of the full cubic group are characterized by one of the irreducible rep- 
resentations of the special cubic group and the parity ±1 which indicates the transformation behaviour 
of the wave functions under the space reflections r — ► —r. These representations will be denoted by 
superscripts Af , A~[ etc. Since the two-particle sector of the Ising model consists of symmetric wave 
functions we will be concerned with even parity representations only. Representations of the cubic 
group naturally arise from representations of the orthogonal group 0(3). This can be illustrated by 
considering the vector space Xi °f an homogeneous polynomials of degree / in r £ IR 3 as the represen- 
tation space of the /-dimensional irreducible representation of 0(3). Through the obvious restriction, 
Xi is also a representation space of the cubic group 0(3, 7L). The decomposition into irreducible repre- 
sentations can be easily found by analyzing the properties of the basic polynomials of Xi under cubic 
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Table 1: The coefficients 8^'. 



k = 0.01 k = 0.05 k = 0.1 





-0.3169 


-0.3247 


-0.3334 


s (2 l 


2.6 10" 5 


5.6 10" 4 


0.0019 




-0.0141 


-0.0157 


-0.0184 


t« 


1.7 10" 5 


4.1 10" 4 


0.0015 


VI 


-0.0074 


-0.0033 


-7.2 10- 


VI 


5.1 10" 5 


2.3 10" 4 


4.1 10" 4 



transformations. For / = 0,2,4 one gets 

O = A+ , 

2 = E+QT+ 

4 = A+QE+QT+® T+ . 

The eigenvalues are parametrized through phase shifts S r as exp (2i<5 r ), where the index T denotes 
the corresponding irreducible representation of the cubic group. At small momentum we make the 
ansatz 

oo 

«r = Eft • (5-29) 

To calculate the coefficients 6^ up to the fifth order we use Bloch's perturbation theory. Since the 
scattering matrix is unitary but not hermitian we introduce 

S r = + = V + Wql , (5.30) 

where 

V = H + S$ q* + S 2 r ql . (5.31) 

From (5.21) and (5.27) it follows that V is diagonal. The non-trivial eigenvalue Vi = 1 + (4/7r)ang* + 
(4/7r)a 2 i^ has multiplicity one, so Vi corresponds to the representation Af . Using Bloch's formula 
(2.15) the higher order contributions can be determined straightforwardly 

45 W W 

cos(2V) + sin(2V) = v 1 + W 1A ql + ^ h " T qt (5.32) 

n = 2 V l~ 1 

To find out if there are more scattering channels at low energies we have to apply perturbation theory 
to the eigenspace corresponding to the infinitely degenerate eigenvalue v 2 = 1. Up to order 0{ql) this 
requires the diagonalization of the reduced matrix 

W)iJ = W i+1J+1 ql + W ^ W ^ ^ . (5.33) 

i — v i 

It is easy to convince oneself that this (44 X 44) matrix is in fact of the order 0{ql). So it can be 
diagonalized numerically. As a result we find one non-trivial eigenvalue of multiplicity two which is 
therefore related to the two-dimensional irreducible representation E + . 

The resulting coefficients S^^ for several values of the inverse coupling k are collected in table 1. 
Since the effective potential is short range one expects the s-wave scattering to be dominant at low 
energy. 
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5.3 Beyond the small momentum regime 

In general, the knowledge about the small momentum behaviour of the scattering quantities is not 
sufficient to answer more subtle questions such as the existence of resonances, for example. Therefore 
the investigation must be extended beyond the small momentum regime. This means that we have 
to find a diagonalization of the scattering matrix independent of the value of the two-particle energy. 
What makes life difficult is the fact that on the lattice there is no rotational symmetry. As a con- 
sequence there is no obvious choice of eigenstates for the scattering operator. Contrary to the small 
momentum expansion the methods will be different for different models and orders of truncation, and 
no general strategy can be presented. 

For the present case, a promising procedure is suggested by the dispersion relation (5.13). Since 
the energy is a simple polynomial of q it seems reasonable to introduce new momentum variables 

q k = 2sm(%) , -2<q k <2. (5.34) 



The plane waves and the (improper) scattering states which are related to & s (q) and $™(q) via the 
coordinate transformation (5.34) are denoted by 3> s (<f) and $™(q), respectively. The normalization 
condition (3.16) translates to 

(Mn *.m = *?m = ^im'-Q) + m'+q)} , (5.35) 

where the normalization function is given by 



Z(q) = y(4- q l)(4-ql)(4-ql) . (5.36) 

We make a further change of coordinates and introduce momentum spherical variables 

q 1 = q s'm(9) cos(^) , q 2 = qs'm(9) sin(^) , q 3 = qcos(9) . (5.37) 

Then eq. (5.35) reads 

A K — Ok 2 ) a 2 — 6"l ~ 

(Mq'), Mq)) = ^ 1 z(E 2 ,n) s(e' 2 -e 2 ) s B (n' - n) , (5.38) 



where 

~ f 1 " V 1 - 2K 2 ml(E 2 - 2m) \ 1/2 



and the normalization function becomes 



(5.39) 



Z(E 2 ,ti) = Z In^^J if " 2 < & [M>*) ^ 2 ' (5.40) 
[ otherwise . 

The next step is to expand the free states 3> s (<f) and the scattering states $™(q) into spherical 
harmonics 



$ s(9) = E E Y 2 ^)\2lm;E 2 ) , (5.41) 

;=0 m = -2f 

oo 21 

^(9) = E E Y£(Si)\2lm;E 2 m) . (5.42) 



;=0 m = -2f 

The new basis vectors are not orthonormal. Instead, from (5.38) it follows that 

(2l'm';E' 2 \ 21 m; E 2 ) = (2l'm';E' 2 in | 21m;E 2 in) = k Z{E 2 ) 2V m , t2lm S(E' 2 -E 2 ) (5.43) 
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where we have introduced the normalization matrix 



z(E 2 ) 2Vm ,, 2lm = 4 2K(g — ^ / dnz(E 2 ,n)Y$(n)Y£(ny . 



(5.44) 



Starting from eq. (3.11) and using the above formulae it is an easy task to expand the scattering 
matrix 

oo 21' oo 21 

S(q',q) = S(E' 2 -E 2 )J2 £ £ £ ^'(fi')* S(E 2 ) 2l . m ., 2lm 



Y 2 7(ty , 



(5.45) 



('=0 m'=-2V i=0 m = 2/ 



with 
and 



2 J2f m',21m 



kZ(E- 



2)21' m' ,21 m 



2vriK T{E 



2 2V m' ,21m 



T(E 



2 21' m'2U 



dfi' /dftY™(iy)T(q',q)Y 2 7(ft)* . 



(5.46) 



(5.47) 



In order to calculate the phase shifts it remains to diagonalize (5.46), i.e. we have to solve the gener- 
alized eigenvalue problem 

Sv = sZv . (5.48) 

Since the normalization matrix (5.44) is invertible, the solution of the eigenvalue equation (5.48) for 
the matrices S and Z is equivalent to the diagonalization of the product 



S'(E 2 ) = [KZiE,]]- 1 S(E 2 ) = 1 - 27TiZ- 1 (E 2 )T(E 2 



(5.49) 



Thanks to the well-suited momentum variables the determination of the transition matrix (5.47) is 
comparatively simple. Since the momentum dependence of the Green's function is via q only, it is just 
the first part of (3.50) that varies with the angular variables We execute the transformation 

to the momentum spherical variables and express the result in terms of spherical harmonics 



J2 cos (lk)cos(q' k ) 



k = l 



^(g 2 -6) 2 + r|^ (^2°(^) 



+t y b {y 2 \Q) + Y 2 - 2 (f2)} {Y 2 2 (fV) + Kf 2 (fV)} . (5.50) 

Using the orthogonality relations the on-shell transition matrix can now be determined straightfor- 
wardly. One obtains 



T(E 



2 21' m' ,21m 



2vr 2 



:{26 l i 1 S m > 6u 6 m0 + <*)/<! Sn (S m > 2 — S m >_ 2 )(S m2 — 6 m _ 2 )} , 



where 



60vr : 



/3(g) = -G(g;0)- 1 + ^ 



(5.51) 



q -6 + 6- 



G(q; 0) 

Finally we substite (5.51) into eq. (5.49) and obtain the scattering matrix 



S'(E 2 



1 



l - { mS [(E 2 )-t^S' 2 (E 2 )} , 



(5.52) 



where 



S[{E 2 ) 2 ii l i2 l m 
> 2 ( 



S',,{E 2 ) 2 ii ; 2 f m 



Z 1 (E 2 ) 2 ji m i oo Sj S m0 

2Z 1 (E 2 ) 2l i m i 20 Sn S m0 

+ [Z 1 (E 2 ) 2l i m i 22 + Z 1 (E 2 ) 2l i m i )2 _ 2 ] [SnS m2 



m-2\ 
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In order to calculate the inverse of the normalization matrix (5.44) one has to distinguish between 
two cases. If q < 2 the integration in (5.44) is over the whole surface of the sphere. Using the 
orthogonality of the spherical harmonics it is straightforward to show that 

z- 1 (E 2 ) 2llmlt2lm = A o * ^ f do z-\E 2 ,h)Ytf \h)Y£(hy . (5.53) 

4 — ZK(q — 0) J 

On the other hand, if q > 2, the multiplication of (5.44) and (5.53) does not lead to the identity 
matrix since the integration is restricted corresponding to (5.40) and the orthogonality relation does 
not apply. In this case we have to calculate Z using eq. (5.44) and invert it afterwards. Because the 
matrices are infinite this procedure requires a truncation of the basis {|2/m; E 2 )} keeping only a finite 
number of elements. For simplicity and since I am not so much interested in lattice artefacts I restrict 
attention to the region q < 2 in the following. The inverse normalization function can be expanded in 
powers of the momentum square q 2 

oo 

Z- 1 (E 2 ,£l) = J2U&)(q 2 y ■ (5.54) 

The coefficients are trigonometric functions of the angular coordinates (9, (f>) and can be ex- 

pressed through spherical harmonics. So the r.h.s. of eq. (5.53) is a sum of integrals over products of 
three spherical harmonics. From the representation theory of the rotation group we know that such 
integrals are related to Wigner's 3j-symbols through (cf. ref. [17], for example) 




/ 47T 

Using the algebraic manipulation programme MAPLE [23] the series expansion (5.54) and the subse- 
quent decomposition of the coefficients Cn(^) into spherical harmonics can be automated completely. 
The non- vanishing 3j-symbols are implemented through Wigner's formula [24] 



h h h \ ( _ 1)h - h - m3 (h + li-l2)Kh-li + h)Kh + l2-h)Kh-m 3 )l(l 3 + m 3 y. 
[ m x m 2 m 3 ) 1 > y (/ 1 + / 2 + / 3 + l)!(/ 1 - mi )!(/ 1 + mi )!(/ 2 - m2 )!(/ 2 + m2 )! 

^ (-l) fc+ ' 2+m2 (/ 2 + / 3 + m 1 -fc)!(/ 1 -m 1 + fc)! 

x tt, kKis-h+h-ky.ih-m^ky.ih-h+m^ky. ' 

where ki = max{0, l\ — 1 2 — m 3 } and k 2 = min{/ 3 — m 3 , l 3 +l 2 — h} . Obviously, the more the momentum 
absolute value approaches q = 2, where Z~ 1 (E 2 ,£l) is singular, the more terms are needed in the 
expansion (5.54). Here we decided to go as far as q = 1.8 which already requires an order of truncation 
of 20 to guarantee a relative accuracy of 1 percent. The CPU time needed was a few hours on an HP 
9000/735 RISC workstation. 

Since the correction to the identity in (5.52) is zero except for the first six columns the diago- 
nalization of the infinite scattering matrix S'{E 2 ) is reduced to the diagonalization of the upper left 
(6 X 6) submatrix of (5.52). As for the small momentum regime we only find two non-trivial eigenval- 
ues of multiplicity one and two corresponding to the scattering phases 6 A + and S E+ . The results for 
k = 1/100 and k = 1/10 as a function of the momentum 

q* = \Jm if (E 2 - 2m) = q y^l - -K 2 m^q 2 . (5.55) 

are plotted in fig. 6. Further scattering channels will only appear if the range of the potential increases, 
i.e. if one includes the higher order contributions to the effective two-particle Hamiltonian. 
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Figure 6: The scattering phases as a function of the momentum for k = 1/100 (dashed lines) and 
k = 1/10 (straight lines). The numerical error is less than 1 per cent for each value. 
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6 Concluding remarks 

The purpose of this paper was to establish a technique to obtain strong coupling expansions for 
scattering phase shifts in Hamiltonian lattice field theories. The method was developed by means of 
the (d+ l)-dimensional Ising model. Starting from Bloch's elegant perturbation theory for degenerate 
levels I derived a convergent series representation for the scattering states and the transition matrix. 
To demonstrate the feasibility of the approach the next-to-leading order approximation for the phase 
shifts in d = 1 and d = 3 spatial dimensions was computed. 

The generalization of the present work to higher orders of the inverse coupling k is straightforward 
although the complexity of the effective two-particle Hamiltonian will rapidly increase. A further 
study of the Ising model might be interesting for its own sake, but I think that it is more promising 
to proceed in another direction. As alredady pointed out in the introduction the lattice Hamiltonian 
of pure gauge theory and QCD has the same general structure (1.1) as the model discussed here. Of 
course I am aware that these theories are much more involved than my test case. In particular the 
static two-particle spectrum is more extensive and will lead to more complicated effective Hamiltonians. 
But hopefully, these technical problems do not mean an insurmountable obstacle towards a qualitative 
insight into the low energy dynamics of hadrons. 

Acknowledgements: The author whishes to thank Martin Liischer for valuable and stimulating 
discussions on the subject, his constant interest in this work and a critical reading of the manuscript. 
Financial support by Studienstiftung des deutschen Volkes is gratefully acknowledged. 

Appendix A 

In this appendix condition (2.41) is checked for the normalizable one-particle states (2.40). To show 
that they are well localized we verify (2.41) in the limit \x\ — ► oo to leading and next-to-leading order 
for the simplest local hermitean operators of the model. The generalization to more complicated, 
composite local observables is staightforward. 

Substituting the explicit formulae for |/) and the ground state (2.21), a straightforward calculation 
leads to 



In case that the fundamental operator (Ji(x) is chosen as local observable, the leading and next- 
to-leading order contributions to both (A.l) and (A.2) vanish. This is because ai(x)\£l)° £ Si and 
<Ji(x)\y)° G {S U S 2 }, so the matrix elements match eigenstates from different eigenspaces of IH . 

Next we consider <D(x) = 03(2:). We remember <7 3 (a;)|fi) = as the defining relation for the 
static vacuum and recognize that <7 2 (a;)|fi) = i|a;) . Using the commutation relations of the Pauli 




(A.l) 



for the l.h.s. and 




(A.2) 



for the r.h.s. 
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matrices one obtains 

<r 3 (x)\y)° = a 3 (x)ai(y)\n)° = {1 - 2S(x-y)}\y)° . (A.3) 

Therefore the next-to-leading order contribution to (A.l) as well as (A. 2) is zero. It remains to 
calculate the leading order of (A.l) 

(f\a 3 (x)\f) = (2vr) d (||/|| 2 - 2\f(x)f) + O(k') . (A.4) 

Since f(x) vanishes in the limit \x\ — ► oo, the localization condition (2.41) is satisfied. 
Finally we want to discuss a less simple situation and choose 

d 

©(a:) = ^ ai(x)ai(x-\-k) . 

k = l 

Here the leading order contribution to (A. 2) is zero because 0(a;)|fi) £ £ 2 , and the sum over the first 
order matrix elements yields 2d. The evaluation of (A.l) is also straightforward and results in 

</|0(*)l/> = \ ^(2vr) d ||/|| 2 + (2vr) d g [r (x; k) + jRi(x;k)} + O(^) , (A.5) 

where the local expressions 

R (x;k) = f*(x + k)f(x) + c.c. 

d 

R^k) = J2{f( x ) [f*(x+k+i) + r(x+k-i)]+ f( x +k) [f(x+i) + f(x-i)]} 
1=1 

-f(x)f*(x)- f( x + k)f(x + k) + c.c. 
vanish in the limit \x\ — ► oo. 



Appendix B 

We show that the in-going scattering states (4.17) of the (1 + l)-dimensional model satisfy condition 
(4.18) for the fundamental operators <D(x) = (Ji(x) and <D(x) = 03(2:). 

The r.h.s. of (4.18) is straightforwardly evaluated using the results of appendix A where the local- 
ization property of the model's one-particle states was established. So for <D(x) = (Ji(x) the leading 
and next-to-leading order contributions vanish. From the explicit form of (4.17) it is easy to see that 
this is also true for the l.h.s. since <Ji(x)\y, z)° £ {Si U £3}, i.e. all matrix elements that occur match 
states from different satic eigenspaces. 

In case of <D(x) = 03(2:) one finds 

r.h.s. = 47r 2 ||/ 1 || 2 ||/ 2 || 2 - Stt 2 { | ^ ( ^ ) | 2 1 1 y 2 1 1 2 + | y 2 ( ^ ) | 2 1 1 ^ 1 1 2 } + O(k') . (B.l) 
For the calculation of the l.h.s. we remark that 

a 3 (x)\y,z)° = {l-2S(x-y)-2S(x-z)}\y,z)° . 
Therefore the next-to-leading order contributions vanish and the result is 

l.h.s. = -L d P ;7 2 (K)*... clp s/ / 1 (^)(^e 1 ^-^)^ + ^sin(|g||i/-^|)sm(|g'|b-z|) 

Z7T J ~ 7r J ~ 7r {y,z 

-4^e^ p - p '^ +y hm(\q\\x-y\)sm(\q'\\x-y\)\ + 0{k 2 ) . (B.2) 
y ) 
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We fix the wave functions fi(p x ) and f 2 (p y ) such that they are supported around momentum values 
of different sign, say p x > and p y < 0. Since 

v d ^ = 2K sin(p) + k 2 sin(2p) + C(k 3 ) . 

the corresponding group velocities are of different sign too, and their ranges do not overlap. It follows 
that q > and q' > in (B.2). So we can forget about the absolute values and write 

e KP-P')(y+*) s - m [ q ( y - z ^ sin[g'(2/-z)] = ± {e'^-^+^-P^ 

_l_ e i[2/(ps-pi)+^(p M -p'j,)] _ e i [y(py-p'J+ z (p^-Py)] _ e i [y(p^-Py)+ z (py-p'J]^ _ (B.3) 

Now the summation can be executed trivially. The momentum ^-functions that result from the last 
two terms of (B.3) do not coincide with the supports of fi and f 2 though. Therefore they do not 
contribute to the subsequent p-integration. Carrying out these steps we find that the final result is 
equal to (B.l). 

In a similar way an explicit check of condition (4.18) in case of ®(x) = ai(x)ai(x-\-l) is possible, 
but requires considerably more elaborate calculations. 

Appendix C 

We discuss the calculation of the three-dimensional Green's function. Our starting point is the integral 
representation 



d 3 p exp(ipr) 



I 3 e 2 {q) - e 2 (p)+ip 



(C.l) 



The next-to-leading order approximation of the two-particle energy is a polynomial in q 

e 2 (q) = -12(1 + k) + 2(1 + 3«)g 2 - ^q 4 . 
For k < k c the denominator of (C.l) factorizes as follows 

e 2 (<7) - e 2 (p) = -| (g 2 " f + v) (f + P 2 - 12 - ^ 
As a consequence we can decompose the Green's function in two integrals of the same structure 

2 - K(q - 6) 



where 



d 3 p exp(ipr) 



\3 -2 ~> . ■ 1 

) q — p l + \p 
d 3 p exp(ipr) 



) 3 q + p 2 - 12 - 4/k 



(C.3) 



(C.4) 



The first step towards the calculation of these three-dimensional integrals is to transform them 
into one- dimensional ones. To begin with we write 

i f°° ( 1 { 2 \] 3 yT 

Gi(q;r) = -- J dtexp ^-t (q 2 - 6) j e~ tp Y[ J — exp(ip k r k ) exp(ii cosp k ) . (C.5) 
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The reader may verify that (C.5) is indeed equal to (C.3) by simply executing the ^-integration. 
Ref. [20, sect. 8.411] tells us that the integrals over the momentum components represent Bessel 
functions of the first kind 

Ci(g;r) = --y o diexpj-^-ejjn^W = I dtF.it). (C.6) 

The ^-integration and the limit p could be interchanged here because the integrand vanishes as 
t~ 3 / 2 for t oo such that the integral (C.6) is well defined. Obviously such a heat kernel representation 
can also be deduced for (C.4) 

C MT,r) = -J o die X p|--^g 2 -6--J|ni r,t ^W = ^ dtF 2 (t) . (C.7) 

The Bessel functions have asymptotic expansions for large arguments [20, sect. 8.451]. This leads to 
a corresponding expansion for the above functions Fi(t) and F 2 (t) Fj(t),j =1,2 

3 oo 

^W = (-i)^EE^W e " bjr " 3/2 > i = i,2, (c.8) 



where 



9 

fi = v = 



Zl = 1(12 - 4/x - f ) , ^ 2 = ^(g 2 - 4/x - 4/«) . (C.9) 

The b^ yir) are analytic functions of the components r fc and the series is rapidly convergent for large 
values of /. 

Next we decompose the integrals (C.6) and (C.7) in three parts 

rtn yoo yoo 

(?,■(«/; 7.)=/ dtFj(t)+ dtFj f (t)+ dtiF^-F^t)} . (CIO) 
Jo Jt N Jt N 

Here F^(t) denotes the asymptotic expansion (C.8) truncated after v = N and t N is some positive 
number. If we substitute (C.8) the integrand of the third contribution to (C.10) is just the rest of the 
asymptotic expansion. One can estimate 

\ Fl (t)- if (0| <\\t\- N -^iZ\b», N+1 (r)\, 

1 /i = 



hence 

/>oo 

dtiF^-Ffit)} 



<^N N - 3/2 i2K N+1 (r)\. (C.ll) 

1 /i = 



We demand that (C.ll) is smaller than some error e and choose some order of truncation N. For 
any r eq. (C.ll) then fixes the value of t N and one obtains an approximation of the Green's function 
where the error is perfectly under control 

r-t N /.OO 

G J (q;r)= dtF^t) + / dtFf(i) + O(s) . (C.12) 
Jo Jt N 

To guarantee an accurracy of 10 digits I chose e = 10~ 12 . The order of truncation was taken to be 
N = 15 which implies t N = 12 for r = and t N = 11 for r = 1. 

In the following we will first describe how the representation (C.12) is used to calculate the Green's 
function for any given values (q; r) and then discuss the series expansion of (C.12) in powers of q. 
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Making use of [20, sect. 3.381] the infinite integral contributing to (C.12) can be expressed as a 
sum of incomplete Gamma functions 

3 N ( ^<^i/2)4 +1/2 l -1 if Zj = 

2^^""'— 1 ,+1/2- 



*« %=o*=o [ v 2 V(-v - 1/2, Zjt N ) otherwise. 



Fortunately the Gamma functions are available in the algebraic manipulation programme MAPLE, 
so we are finished with this part of (C.12). Concerning the calculation of the finite integral we have 
to take into account that F 2 (t) strongly oscillates with a frequency proportional to 1/k. From several 
tests we found that in this case the accurracy of the numerical integration routine implemented in 
MAPLE is not sufficient unless one subdivides the integration intervall in such a way that each part 
roughly contains one period of the function F 2 (t). This would mean to solve approximately 10/k 
integrals numerically taking a disproportionately large amount of computer time. Therefore, instead 
of a numerical integration routine we use the series representation [20, sect. 8.440] for the Bessel 
functions 

3 / -i \ r 2 +r 2 +r 3 oo / f \ 2n 

nw'>=( 5 ) x>G) (c.i4) 

with some coefficients £„. Since this series is convergent it can be approximated by truncation. If we 
demand that the error is less than e the order of truncation is about 100 which is certainly no problem 
if one has an algebraic manipulation programme at hand. The remaining integrations can be done 
analytically, so there is no further loss of accurracy. 

In case one is only interested in the small momentum behaviour of the Green's function it suffices 
to expand Gj(q; r) in powers of q. Starting from the representation (C.12) such a series can be easily 
derived. Concerning the finite integrals we simply expand the exponentials (cf. eqs. (C.6) and (C.7)) 
and do the ^-integration as described above. To obtain a small momentum expansion for the infinite 
integrals one has to distinguish between two cases. If Zj(q = 0) / we can make use of the lower 
line of eq. (C.13). According to [20, sect. 8.354] the incomplete Gamma functions have the series 
representation 



T(-u -1/2, zjtff) = V{-v -l/2,Xj) 

n!r(i/ + 3/2) ~ / 



,- 3 / 2 f ^-l)"r(i/ + n + 3/2) ^ y] 



n=0 • V 1 / / f=n+l 

where 

Vi = (-!)' g tN ^ ' x i = z i ~ y i ■ ( C - 16 ) 
However, for fj, = 3 it follows that z x = —(i/2)q 2 leading to the contributions 

"oI>» / d*exp{(i/2)^ 2 }r'- 3 / 2 

Z ,,-n Jt N K J 



(C.17) 



Iv + l 



These integrals cannot be represented through incomplete Gamma functions if q = (cf. the upper line 
of eq. C.13). Therefore the expansion around q = must be done differently. By partial integration 
one deduces a recurrence relation 

I v+1 = (u+ 1/2)" 1 (exp {(i/2)^ 2 } r "- 1 ' 2 + X -f /,) . (C.18) 

The first contribution can be expanded trivially whereas the integral is iterated again. What remains 
is the integral 

I = [°° dtexp{(\/2)tq 2 } r 1/2 . (C.19) 
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Table 2: Results for G„(0). 





V 


K = 0.01 


k = 0.05 


K = 0.1 







-0.1238.. 


-0.1149 


-0.1055.. 




1 


-0.0386.. i 


-0.0345.. i 


-0.0306.. i 




2 


-0.0065.. 


-0.0078.. 


-0.0089.. 




3 


-0.0050.. i 


-0.0050.. i 


-0.0050.. i 




4 


0.0075 + 0.0050.. i 


0.0066 + 0.0045.. i 


0.0056.. + 0.0040.. i 




5 


-0.00068.. i 


-0.00070. .i 


-0.00071.. i 






Table 3 


: Results for G V {1). 




V 




K = 0.01 


k = 0.05 


K = 0.1 







-0.0417.. 


-0.0373.. 


-0.0329.. 


1 




-0.0386.. i 


-0.0345.. i 


-0.0306.. i 


2 




0.0143.. 


0.0122.. 


0.0102.. 


3 




0.0014.. i 


0.00068.. i 


0.00009.. i 


4 




-0.0065 + 0.0024.. i 


-0.0051 + 0.0022.. i 


-0.0039.. + 0.0019.. i 


5 




0.00014.. i 


0.00014..i 


0.00011.. i 



If we shift the lower limit of integration to zero we can use [20, sect. 3.691] and obtain 

/„ = (1 - i) - f N diexp \a/2)tq 2 } r 1 ' 2 . (C.20) 
q Jo 1 } 

Collecting the results one ends up with the desired small momentum expansion for the Green's 
function 

oo 

G(q- 1 r) = Y,G v (r)q • (C.21) 

The first five coefficients for r = 0, r = 1 and some values of the inverse coupling k are listed in table 
2 and 3. 



34 



References 

[1] M. Liischer, On a relation between finite size effects and elastic scattering processes, in: Progress 
in gauge field theory, (Cargese 1983), eds. G. 't Hooft et al. (Plenum, New York, 1984) 451; 
M. Liischer, Comm. Math. Phys. 104 (1986) 177 ; 
M. Liischer, Comm. Math. Phys. 105 (1986) 153 ; 
M. Liischer, Nucl. Phys. B354 (1991) 531 

[2] I. Montvay and P. Weisz, Nucl. Phys. B290 (1987) 327 ; 

Ch. Frick, K. Jansen, J. Jersak, I. Montvay, P. Seuferling and G. Miinster, Nucl. Phys. B331 
(1990) 515 

[3] M. Liischer and U. Wolff, Nucl. Phys. B339 (1990) 222 
[4] J. Nishimura, Phys. Lett. B294 (1992) 375 

[5] H.R. Fiebig and R.M. Woloshyn, Meson-meson scattering in lattice QED 2+1 , in : Lattice 92, Eds. 
J. Smit and P. van Baal, Nucl. Phys. B(Proc. Suppl.) 30 (1993) 883 ; 

H.R. Fiebig, R.M. Woloshyn and A. Dominguez, Meson-meson scattering in (2 + l)-dimensional 
lattice QED, Florida Int. Univ. preprint, FIU-PHY-92- 10-23 

[6] M. Liischer, Nucl. Phys. B364 (1991) 237 

[7] C.R. Gattringer and C.B. Lang, Phys. Lett. B274 (1992) 95 ; 
C.R. Gattringer and C.B. Lang, Nucl. Phys. B391 (1993) 463 ; 

C.R. Gattringer, I. Hip and C.B. Lang, Scattering in a simple 2-d lattice model, in : Lattice 92, 
Eds. J. Smit and P. van Baal, Nucl. Phys. B(Proc. Suppl.) 30 (1993) 875 

[8] F. Zimmermann, J. Westphalen, M. Gockeler and H.A. Kastrup, Unstable particles in finite 
volume: The broken phase of the 4-d 0(4) non-linear cr-model, in: Lattice 92, Eds. J. Smit and 
P. van Baal, Nucl. Phys. B(Proc. Suppl.) 30 (1993) 879 

[9] S.R. Sharpe, R. Gupta and G.W. Kilcup, Nucl. Phys. B383 (1992) 309 ; 
R. Gupta, A. Patel and S.R. Sharpe, Phys. Rev. D48 (1993) 388 

[10] Y. Kuramashi, M. Fukugita, H. Mino, M. Okawa and A. Ukawa, Phys. Rev. Lett. 71 (1993) 2387; 
Y. Kuramashi, M. Fukugita, H. Mino, M. Okawa and A. Ukawa, Lattice QCD calculation of 
hardron scattering lengths, KEK preprint, KEK-CP-010 (1993) 

[11] V. Alessandrini and A. Krzywicki, Phys. Lett. 93B (1980) 173 ; 
V. Alessandrini and A. Krzywicki, Phys. Rev. D24 (1981) 3237 

[12] V. Hakim, Nucl. Phys. B180 (1981) 417 

[13] G. Miinster, Nucl. Phys. B249 (1985) 659; 

G. Miinster, The strength of glueball interactions, in: Ahrenshoop Symp. 1985, 248 

[14] E. Fradkin and L. Susskind, Phys. Rev. D17 (1978) 2637 

[15] C. Bloch, Nucl. Phys. 6 (1958) 329 

[16] M. Liischer, Selected topics in lattice field theory, in: Fields, Strings and Critical Phenomena 
(Les Houches 1988, Session XLIX), eds. E. Brezin, J. Zinn- Justin (Elsevier, Amsterdam, 1990) 
451 

[17] A. Messiah, Quantum mechanics, vol. II (North-Holland, Amsterdam, 1962) 



35 



[18] M. Reed and B. Simon, Methods of modern mathematical physics, vol. I - III (Academic Press, 
New York, 1972) 

[19] P. Pfeuty, Ann. Phys. 57 (1970) 79 

[20] I.S. Gradshteyn and I.M. Ryzhik, Table of integrals, series, and products (Academic Press, New 
York, 1980) 

[21] J.B. Kogut, Rev. Mod. Phys. 51 (1979) 659 

[22] M. Sato, T. Miwa and M. Jimbo, Proc. Japan Acad. 53, Ser. A (1977) 6; 
B. Berg, M. Karowski and P. Weisz, Phys. Rev. D19 (1979) 2477 

[23] B.W. Char et al., MAPLE V library reference manual and language reference manual (Springer, 
New York, 1991) 

[24] E.P. Wigner, Group theory (Academic Press, New York, 1959) 



36 



